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Abstract: This article translates some recent results on quasilocal horizons into the language 
of (3 + 1) general relativity so as to make them more useful to numerical relativists. In 
particular quantities are described which characterize how quickly an apparent horizon is 
evolving and how close it is to either equilibrium or extremality. 

1. Introduction 

In (3 + 1) -dimensional numerical relativity an apparent horizon in a three-dimensional spacelike 
slice St is the outermost smooth two-surface whose whose outward null expansion vanishes [1,2]. 
Then an apparent horizon world-tube is a three-surface H foliated by two-surfaces St C with 
vanishing outward null expansion. Such surfaces have been an active topic of research in mathematical 
relativity for the last 10-15 years (see [3-5] for reviews or the reference section of [6] for more recent 
papers) and results from that work are now beginning to be applied to numerical relativity [7, 8]. 

This short paper continues this tradition by translating results from recent papers into the lan- 
guage of (3 + 1) -dimensional relativity. We first consider the recent work on slowly evolving hori- 
zons [6, 9, 10]. These papers found conditions under which a horizon can be considered to be in a 
quasi-equilibrium state. Intuitively, slowly evolving horizons are "almost" isolated [11] and so "nearly" 
null. Here we shall focus on an expansion parameter that arose from that work. This parameter is nec- 
essarily small for slowly evolving horizons, but even if it isn't small it can be used to invariantly 
characterize the rate of expansion of horizons. It is to this end that we discuss it here, with a particular 
eye to tracking how quickly a black hole settles down to equilibrium after a merger or formation. 

Second we consider a parameter that characterizes how close an apparent horizon is to being ex- 
tremal [12]. Extremality in this case is a generalization of the well-known bound on the maximum 
angular momentum of a Kerr black hole relative to its mass (or equivalently surface area). For apparent 
horizons this turns out to be related to the expansion parameter and places a bound on the maximum an- 
gular momentum relative to the intrinsic geometry of the horizon. The immediate use for this parameter 
would be to study the proximity to extremality of post-merger black holes. 

2. Horizons in (3+1)-dimensional gravity 

Let (Af, g^ft, Vq) be a four-dimensional spacetime which is foliated by spacelike three-surfaces 
{(St, /lij , Di) ]ti < t < where hij is the induced metric and Di is the compatible covariant 
derivative. Note the indices; throughout this paper we will use early-alphabet latin indices for four- 
dimensional quantities and mid-alphabet latin indices for those defined in the three-slices. We use e to 
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denote the pull-back/push-forward operator between the various spaces. 

Next let be a compatible time-evolution vector field generating a flow which maps the St into 
each other. That is £ = 1- In the usual way this may be broken up as 

T" Nu" + V , (1) 

where N is the lapse function, = e°T^' is the spacelike shift vector field for some G TEt, and 
u° is the future-pointing timelike normal vector field to the S^. 

Finally we choose our sign convention so that the extrinsic curvature Kij of the St in M is 

Kij = j^{hij - C vhij) ^ e'^e^\IaUb , (2) 

where the dot indicates a Lie derivative with respect to T". 

Now, let St be a closed two-surface in a slice St and write its outward-pointing unit normal as P. 
Then future-outward and future-inward pointing null normals to St in M can be written as 

r = f{u^ + n and - ^ - n , (3) 

where s"' = e°s' and / is an arbitrary positive function: we have adapted the standard cross-normalization 
^ • n = — 1 so there is the only free parameter in their definition. For / = 1 we denote this pair of null 
normals as (^°, n°) as shown in Fig. 1. 

We write the induced two-metric on St as qab (using capital Latin indices for such two-dimensional 
quantities) and note that 

q"'' = e^e'^Bq^^ = .g"^ + + i'^n'' and q'^ = e^eJ^q"^^ = - s's^ . (4) 
Then the outward and inward null expansions are, respectively, 

0(1) = r'^ah = f{q''K,j + If'Dis,) and 0(„) ^ r'^^anb = ^(^^'^y - g'^ As,) • (5) 

As noted in the introduction, on a numerical apparent horizon — and in practical simulations 
where there are known to be black holes in the initial data, this will usually be a sufficient condition to 
identify a black hole boundary. A key word in the previous sentence is "usually" and we will return to 
the potential complications in section 4. 
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For now, however, we will assume that an apparent horizon finding algorithm [2] has been used 
to identify 6'(£) = surfaces on each Et over some range of t and further that over that range their 
union is a smooth three-surface H. Then we can always find a horizon evolution vector field V that is 
normal to the St and maps these surfaces into each other Since the St C we can write 

V" = A^(m" + i)_LS°) , (6) 

where N is again the lapse and v± is the the "velocity" of the horizon relative to the foliation. For a 
null horizon parallel to £"-,v± = 1 while for a spacelike and expanding horizon v± > 1. 

3. Expansion 

We now consider the first of the physical characteristics: the rate of expansion of the horizon. 
Relative to the foliation, the rate of change of the area element ^/g is 

where we have made use of the fact that 6'(£) = <^ (p-' Kij = —q^^DiSj (which could also be used to 
rewrite the last line in terms of q^^ Kij). This expansion is manifestly dependent on the lapse N and so 
the coordinate t labelling of the hypersurfaces E^. However, if H is not null there is an obvious remedy 
to this problem : calculate the rate of expansion with respect to the unit V"^ = V°/||V|| instead of V 
itself. Then the expansion relative to V is 



1 



1 



q''D,s, , (7) 



where a = sgn(w± — 1). Note that Bfy-^ goes to zero as v±^ 1 even though V" itself diverges. 

To better understand ^'(■i>), we return to four-dimensions and the formalism developed for slowly 
evolving horizons [6, 9]. There, V° was written in terms of the null normals and a parameter C: 

V° = r - Cn" . (8) 

Then the foliation of the horizon fixes both a scaling of the null vectors and C. Comparing with (6) it 
is straightforward to see that in the (3 + l)-formalism these take the form: 

f = N{1 + vj_)/2 and C = N'^{vI - 1)/2 . (9) 

The value of C fixes the signature of the horizon: C>0<^|tij^|>l<=>i7is spacelike, C — ■i^ 
\v±\ = 1 -i^ H is null, and C < \v± \ < 1 'i^ H is timelike. We can rewrite 

= -aVT^^(n) , (10) 

where a = sgn(C). In this form 0^yj (or actually its square) is the most important of several quan- 
tities used to decide whether or not a horizon is "almost" isolated and so in quasi-equilibrium with 
its surroundings. Among other properties such horizons will obey approximate zeroth law and first 
laws of black hole mechanics with all fluxes being calculated relative to the normals (as if they 
where truly null). Some intuition for how dramatically a horizon can evolve and still be counted as in 
a quasi-equilibrium state can be gained from the examples of [10]. 

To be mathematically certain about such a classification one needs to track 6^^^ and the other 
quantities point-by-point on the horizons. However to get a quick idea of the rate of expansion (and 
proximity to equilibrium) it is convenient to have a single, dimensionless, number to calculate. Such a 
quantity can be obtained by integrating the square of 6^^^ over St'. 

-1 



area 3 ./^" V . 1 - r („) ^ 



VJ_ + 1 



{q'W^s.y. (11) 



In spherical symmetry this is proportional to {dRn / dsY , where Rh = y/a/An is the areal radius of 
St and s is the arclength "up" the horizon along a flow fine of V". 
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Fig. 2. A horizon H with 6'(£) = and < evolves between equihbrium states. Here, null horizons have 
slopes of ±45°, spacelike horizons have slopes from —45° to 45°, and timelike membranes have the steeper slopes 
outside of that range. Thus before A, H is in equilibrium and null (parallel to i!"). From ^ — > C it is expanding and 
spacelike. At C it is momentarily null (parallel to n") and earea diverges. From C ^ D, is a timelike membrane 
before becoming null again at V. Finally from T) ^ T A'l?, again a dynamical horizon before becoming isolated 
and quiescent after T. Between C and T), H is super-extremal while at those points it is extremal. 



4. Apparent horizons as black hole boundaries 

We now return to the "usually" that followed equation (5) and consider when we can definitively 
associate a black hole with an apparent horizon. Following Penrose [14] we take the existence of fully 
trapped surfaces (0(£) < 0, 0(„) < 0) as being the key characteristic of black holes. Together with the 
null energy condition (and some more technical assumptions) the existence of trapped surfaces both 
implies a spacetime singularity somewhere in their causal future and, if the spacetime is asymptotically 
flat, an enveloping event horizon [13]. Then, taking the set of all points that lie on trapped surfaces as 
the interior of a black hole, one would intuitively expect the boundary of that region to be foliated 
by two-surfaces St with 9(^^) = and 0(„) < 0, and further expect that slight inwards deformations 
of the St should turn them into fully trapped surfaces. These ideas are the foundation for both the 
mathematical relativity definition of an apparent horizon (the closure of the boundary of the set of all 
points lying on trapped surfaces in a given Ej) [13] as well as Hayward's trapping horizons [15]. 

Though in reality things turn out to be more complicated (see for example the discussion of "wild" 
trapped surfaces in [16] or the counter examples in [17]) the motivation remains valuable and = 
"boundaries" with trapped surfaces just inside are, at the very least, associated with black holes and 
their boundaries. In preparation for discussing their properties we quickly review deformations (for 
more details see [6]). A normal deformation to a surface St is characterized by a deforming vector field 
X° G T±St (the normal bundle to St). Then taking any smooth extension of X" into a neighbourhood 
of St one can construct the corresponding flow and use this to deform St ; schematically St ^ St + eX. 
The deformation operator 5x then calculates the rate of change of a quantity during such an evolution. 
This operator is linear in the sense that Sx+y = 5x + Sy but in general S fx 7^ fSx- 

Thus, for an apparent horizon, the assumption that there be trapped surfaces "just inside" St cor- 
responds to there existing an inward pointing vector field X"- E TY^t such that 6xS(e) < 0. If such 
a two-surface exists on one slice Y,t„ then this two-dimensional apparent horizon can necessarily be 
extended into a three-dimensional, apparent horizon world-tube H over some finite range of t [18]. If 
the null energy condition also holds, then H is necessarily spacelike or null (that is C > 0) and non- 
decreasing in area if 6'(„) < [15, 18, 19] — this last fact follows directly from SvVq = —y/qC9(^n)- 



©2008 NRC Canada 



I. Booth 



Null regions are isolated horizons [11] and may be regarded as black hole equilibrium states while the 
spacelike and expanding regions are dynamical horizons [19]. On any given S't , V° will be either space- 
like everywhere or null everywhere [18]. That is, transitions between isolated and dynamical horizons 
must happen "all at once": no individual St can be partly isolated and partly dynamical. 

It is possible for a horizon satisfying 6'(£) = 0, 6'(„) < 0, and 6nd{e) < (note the use of the 
inwards null direction for cases where there isn't a favoured inward spacelike direction) to evolve into 
a timelike membrane on which 6nd{e} > [2,20,21]. Such a situation is shown in Fig. 2. Physically, 
evolutions of this type correspond to situations where new horizons form outside of existing ones. The 
cited examples were for spherically symmetric spacetimes but it is widely speculated that something 
similar happens during at least some of the apparent horizon "jumps" seen in numerical relativity. For 
example, if one was only tracking the outermost Of^^-j =0 surface in Fig. 2, it would jump at t = t^. 

It is useful to compare C and v± in this figure. While the sign of C faithfully tracks the signature 
of H, v± is coordinate dependent and diverges at B and £ even though nothing physically significant 
happens there. For v±, the physically significant values are ±1 not zero or infinity. 

5. Evolutions and extremality for dynamical horizons 

In the example depicted in Fig. 2, C diverges and H changes signature when Sn6(^g-^ = 0. In 
spherical symmetry Sn0{e) doesn't depend on the scaling of the null vectors, but in general one must 
select a "correct" scaling to observe this. For dynamical horizons this scaling is easy to find. Relative 
to the null vectors of section 3, rescale so that i = i/C and n = Cn or equivalently choose / — 
l/{N{v± — 1)). Then since 6^i) = everywhere on H we have [6, 12] 

Svei = CSiO^i) - SnOf^i) ^ where 5i6i ^ -af^ (T(^i)AB - SnGTabre^ . (12) 

^{l)AB — b'a^'b'^ ah is the shear in the Z° direction (note its simplified form thanks to df^f-^ = 0). One 
of the reasons for the scahng choice is now clear: for general deformations Sfx ^ fSx but is a 
special case where multiplicative factors can be commuted through S. From this result it is easily seen 
that if the null energy condition holds and di9(^i^ ^ then C > bnQ(j^ < (and vice versa). This 
is one proof of the already mentioned result that dynamical apparent horizons with trapped surfaces 
just inside are generically spacelike. Further note that if SnO^g-^ then C necessarily diverges. 

In [12] it has been argued that the most natural way to define extremality for dynamical horizons 
is in terms of SnO^^i-^^ . A sub-extremal horizon is one for which Sn^i < (and so there are trapped 
surfaces just inside) while for an extremal horizon Sndi = 0. This definition follows by extension from 
isolated horizons such as Kerr where it implies the usual restrictions on surface gravity (positive for 
sub-extremal, vanishing for extremal) and, in cases where it is well-defined, angular momentum. Note 
however that this characterization also applies to highly dynamical and distorted horizons where the 
angular momentum and surface gravity may not be so well-defined. 

For a dynamical numerical apparent horizon Sn6/i can be calculated directly as it was in [12] but in 
most cases it is probably easier to combine equations (9) and (12) to obtain: 

SnO^i) = (!!r^) (Hio)\\' + »^GTa,eyl) (O) 

where = u"" + s° and so cr[i^)AB ~ ^A^Bi-^ij + DiSj). Clearly this expression goes to zero if 
v± ^ —1 (H is null and parallel to — n") and becomes positive for timelike membranes. As for the 
expansion we can define a dimensionless parameter that tracks how close a horizon is to extremality. 



'However in near-isolated cases where C ~ and so the suggested reseating £ — £/C may be numerically 
inappropriate, other methods of calculating this quantity may be more useful [12]. 
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We call this the extremality parameter. 

6 = 1 + ^/ dPx^qS^e^,) = 1 - / d^x^q (^1^^) (I I ? + SnGTatrjl) . (14) 
^ ^ St ^ ^ St V ^'-L / 

If the dominant energy condition holds e > (this is most easily seen from the alternative method of 
calculating quantity described in [12]). Thus, for sub-extremal holes < e < 1 while for an extremal 
horizon e = 1. Super-extremal time like membranes will have e > 1. Note however that the intended 
use of this parameter is not to determine the signature of the horizon but rather to quantify how close 
a horizon is to extremality and so a transition to a timelike membrane. If one is only interested in 
signature, v± is sufficient. 

6. Conclusions 

This short paper has defined two dimensionless quantities earea (Eq. 11) and e (Eq. 14) which can 
be used to characterize dynamic black hole evolutions, earea invariantly defines the rate of expansion: 
it is very small for slowly evolving horizons and vanishes only for truly isolated horizons. As such it is 
well suited to tracking the approach of a dynamical horizon to equilibrium. By contrast e is probably 
of greatest interest in highly dynamical situations. As noted < e < 1 for regular sub-extremal 
holes, achieves unity for extremal horizons and becomes greater than one for super-extremal timelike 
membranes. Thus it is well suited to tracking the approach to or aftermath of transitions between 
timelike membranes and dynamical horizons — though it should be kept in mind that earea also has 
something to say about such situations as it diverges whenever e = 1. 
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